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ABSTRACT 

We describe the angular power spectrum of resolved sources at 3.6 pm (L-band) in Spitzer imaging data of the 
GOODS HDF-N, the GOODS CDF-S, and the NDWFS Bootes field in several source magnitude bins. We also 
measure angular power spectra of resolved sources in the Bootes field at K$ and /-bands using ground-based 
IR imaging data. In the three bands, J, K$, and L, we detect the clustering of galaxies on top of the shot- 
noise power spectrum at multipoles between £ ~ 10 2 and 10 5 . The angular power spectra range from the large, 
linear scales to small, non-linear scales of galaxy clustering, and in some magnitude ranges, show departure 
from a power-law clustering spectrum. We consider a halo model to describe clustering measurements and 
to establish the halo occupation number parameters of IR bright galaxies at redshifts around one. The typical 
halo mass scale at which two or more IR galaxies with L-band Vega magnitude between 17 and 19 are found 
in the same halo is between 9 x 10 11 M Q and 7 x 10 12 M Q at the la confidence level; this is consistent with 
the previous halo mass estimates for bright, red galaxies at z ~ 1 . We also extend our clustering results and 
completeness-corrected faint source number counts in GOODS fields to understand the underlying nature of 
unresolved sources responsible for IR background (IRB) anisotropics that were detected in deep Spitzer images. 
While these unresolved fluctuations were measured at sub-arcminute angular scales, if a high-redshift diffuse 
component associated with first galaxies exists in the IRB, then it's clustering properties are best studied with 
shallow, wide-field images that allow a measurement of the clustering spectrum from a few degrees to arcminute 
angular scales. 

Subject headings: cosmology: theory — large scale structure of universe - 
ies 



- diffuse radiation — infrared: galax- 



1. INTRODUCTION 

The intensity of the cosmic near-infrared background (IRB) 
is a measure of the total light emitted by stars and galaxies 
in the Universe. While the absolute background has been 
estimated by space-based experiments, such as the Diffuse 
Infrared Background Experiment (DIRBE; Hauser & Dwek 
2001) and the Infra-Red Telescope in Space (IRTS; Mat- 
sumoto et al. 2005), the total IRB intensity measured still re- 
mains fully unaccounted for by sources: only 13.5 ± 4.2 nW 
m" 2 sr" 1 is resolved to point sources at 1.25 /im (Cambresy et 
al. 2001), while current direct measurements range from 25- 
70 nW m" 2 sr" 1 . At wavelengths greater than 3 /im, with a to- 
tal IRB intensity of 12.4 ± 3.2 nW m" 2 sr" 1 (Wright & Reese 
2000) at 3.6 ^m, the "missing source" problem is less signifi- 
cant with 5.4 nW m" 2 sr" 1 resolved to point sources in Spitzer 
data (Fazio et al. 2004a; Franceschini et al. 2006), while fluc- 
tuation analyses of Spitzer data yield 10.6^ 95 nW m" 2 sr" 1 
(Savage & Oliver 2005). We refer the reader to Kashlinsky 
(2005) for a recent review on the IRB including a summary of 
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past attempts to understand the intensity excess relative to the 
background predicted by resolved source number counts. 

Primordial galaxies at redshifts 8 and higher, especially 
those involving Population III stars, are generally invoked 
to explain the missing IR flux between 1 and 2 pm, 
with most of the intensity associated with redshifted Lyman- 
a emission during reionization (e.g., Santos, Bromm & 
Kamionkowski 2002; Salvaterra & Ferrara 2003; Cooray & 
Yoshida 2004; Fernandez & Komatsu 2006). While models 
of high-redshift Pop III populations can explain the "missing" 
IRB, these models run into several difficulties if such sources 
were to account for all of the missing IR intensity. These in- 
clude the high efficiency required to convert baryons to stars 
in first galaxies (e.g., Madau & Silk 2005) and limits from 
deep IR imaging data that suggest a lack of a large population 
of high-redshift dropouts (Salvaterra & Ferrara 2006). The re- 
cently revised optical depth to reionization (Page et al. 2006), 
with reionization around a redshift of 10 instead of the pre- 
vious estimates of 20, decreases the number density of first 
sources that are required for reionization. In return, the im- 
plied fractional contribution from first sources to the total in- 
tensity of the IRB is lowered. 

Still, one does expect some contribution to the IRB from 
sources that reionized the Universe, though the exact inten- 
sity of the IRB from such sources is yet unknown both the- 
oretically and observationally. As pointed out in Cooray 
et al. (2004; also, Kashlinsky et al. 2004), if a high- 
redshift population contributes significantly to the IRB, then 
these sources are expected to leave a distinct signal in the 
anisotropy fluctuations of the near-IR intensity, when com- 
pared to the anisotropy spectrum associated with low-redshift 
sources. Using results from a fluctuation analysis in deep 
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Spitzer imaging data with resolved point sources removed 
from the image down to a deeper magnitude than previous 
studies on this topic, Kashlinsky et al. (2005) claimed a po- 
tential detection of the clustering signature of high-redshift 
sources at wavelengths of 3.6, 4.5, and 5.8 fim. A previous 
attempt to understand the nature of this excess clustering in 
Spitzer data, when resolved sources are removed down to a 
magnitude level of 22.5 in the L-band, suggested that it could 
be the clustering signature of galaxies at redshifts greater than 
5, with a total contribution to the IRB at the level of ~ 1 to 2 
nW m" 2 sr" 1 in the L-band (Salvaterra et al. 2006); at shorter 
wavelengths, the intensity of this background remains uncer- 
tain as the spectrum of this excess fluctuation component is 
only established with Spitzer for A > 3.6 /im. While this sug- 
gestion was simply based on a model description of the fluc- 
tuation spectrum, whether such a scenario is consistent with 
faint source counts in deep Spitzer images is yet to be consid- 
ered. 

If the 1 to 2 /im IRB is not fully resolved by extragalactic 
sources, then it could be that most of the missing source prob- 
lem is due to deficiencies in accounting for the zodiacal light 
(Dwek, Arendt & Krennrich 2005). This possibility may be 
further consistent with some suggestions that, subject to a pri- 
ori assumptions on their intrinsic spectral shape, the observed 
blazar TeV spectra rule out a significant IRB from extragalac- 
tic sources (e.g., Aharonian et al. 2006). While the point 
source photometry may underestimate the total IR flux in each 
resolved source by missing the outer regions when construct- 
ing source catalogs, it is unlikely that this effect alone can 
explain the factor of 2 to 3 difference between predicted and 
measured IRB intensity given that point source photometry is 
unlikely to miss more than 5% to 10% of the total flux (see 
Totani et al. 2001 for a discussion). Such a small residual, 
however, could produce an excess in fluctuations of the IRB, 
with the angular spectrum tracing that of the bright resolved 
sources. 

To separate various possibilities related to the nature of 
intensity as well as anisotropy excesses, further analyses of 
near-IR background data is required. If the excess background 
is due to missing flux in low-z galaxies, then one would ex- 
pect fluctuations in the background to strongly correlate with 
resolved point sources. The contribution from zodiacal light 
is known only to within an upper limit, but if residual zodia- 
cal light is the reason, then spatial clustering of the IRB will 
peak at degree or more angular scales due to its smoothness. 
Thus, to distinguish between these possibilities, it is necessary 
to further understand spatial clustering at IRB wavelengths. 
Here, we concentrate on the fluctuations in IRB due to re- 
solved sources. We study both the shape and the amplitude 
of the angular power spectrum and use these measurements to 
make a relative comparison between resolved and unresolved 
IRB clustering, with the latter from measurements in Kash- 
linsky et al. (2005). We also connect our clustering measure- 
ments with faint number counts and an underlying halo model 
of the IR galaxy distribution within the large-scale structure to 
study if faint galaxies could be responsible for the IRB fluc- 
tuations seen with Spitzer. 

In this paper, we make clustering measurements at 3.6 fim 
in the Great Observatories Origins Deep Survey (GOODS; 
Dickinson et al. 2003) HDF-N field, the GOODS CDF-S 
field, and the NOAO Deep Wide-Field Survey (NDWFS; Jan- 
nuzi & Dey 1999) Bootes field with the source sample divided 
into several magnitude bins in each of these catalogs. The data 
come from Spitzer IRAC imaging of these fields. We refer the 



reader to Dickinson et al. (in preparation) for IRAC observa- 
tions of the GOODS field and Eisenhardt et al. (2004; also 
see Fazio et al. 2004a) for the Bootes-field IRAC shallow sur- 
vey. In addition to Spitzer images, we also measure clustering 
of resolved sources in the Bootes field in ground-based im- 
ages at K s and /-bands with data taken by the FLAMINGOS 
Extragalactic Survey (FLAMEX; Elston et al. 2006). 

Here we present clustering measurements in terms of an 
angular power spectrum in multipole space. Due to simplic- 
ity in the measurement, galaxy clustering is also measured in 
real space in terms of the angular correlation function or w(9). 
As discussed in Tegmark et al. (2002), Fourier or multipole- 
space measurements should be preferred over the real space 
due to the simple behavior of the covariance matrix and the 
fact that multipole-space measurements can be easily related 
to the underlying three-dimensional power spectrum without 
a complex window function that maps fluctuations in Fourier 
space to that of the angular space. In the case of the real space 
correlation function, clustering at different angular scales are 
also correlated (Eisenstein & Zaldarriaga 2001) and this also 
lead to a complex covariance matrix that must be accounted 
for when model fitting the measurements. 

Despite these considerations, real-space clustering mea- 
surements are easy to make and have been studied in the lit- 
erature at 3.6 fim and other Spitzer IRAC bands (Oliver et al. 
2004; Fang et al. 2004). While model fitting was considered 
in terms of a power-law, a clear departure from a power-law 
model fit is present in angular clustering of 3.6 /im galaxies 
in the Spitzer Wide-Area Infrared Extragalactic (SWIRE) as 
described by Oliver et al. (2004; for example, see their Fig- 
ure 1 for clustering in Lockman and EN1 fields that indicates 
a rise above a power-law at small angular scales). On the 
other hand, clustering measurements in the First Look Survey 
(FLS) at 3.6 /im by Fang et al. (2004) do not show clear ev- 
idences for a departure from the power-law (see, Figure 1 in 
Fang et al. 2004). This could be due to the smaller number of 
galaxies used for clustering measurements, when compared 
to statistics in Oliver et al. (2004). Here, using the Bootes- 
field IRAC shallow survey over 6 deg. 2 , we find a departure 
from the power-law description and model that departure in 
terms of the halo model for galaxy clustering (Cooray & Sheth 
2002). A departure from a power-law description, however, is 
not surprising for galaxy clustering at high redshifts due to 
non-linear clustering at angular scales less than a few arcmin- 
utes (Ouchi et al. 2005; Lee et al. 2005). Such departures 
have been detected for Lyman-break galaxies at z ~ 3 to 5 
(Ouchi et al. 2005; Lee et al. 2005; Cooray & Ouchi 2006), 
as well as at low redshifts with surveys such as SDSS (Ze- 
havi et al. 2005), DEEP2 (Coil et al. 2004), and COMBO-17 
(Phleps et al. 2004). These departures were expected, and 
are described, in terms of analytical approaches such as the 
halo model (Cooray & Sheth 2002) and conditional luminos- 
ity functions (Cooray 2006a; van den Bosch et al. 2006). 

The Paper is organized as following: in the next section, we 
briefly summarize our data and in § 3 we describe statistics 
in our catalogs. In § 4, we describe the procedure we used 
to measure clustering of IR sources and in § 5 we present a 
summary of the analytical model used to describe clustering. 
Based on model fits, we extract quantities related to how IR 
galaxies occupy dark matter halos. In § 6, we discuss our re- 
sults in the context of recent IRB anisotropy measurements 
with Spitzer and future wide-field surveys with Akari (Mat- 
suhara et al. 2006), Cosmic Infrared Background Explorer 
(CIBER; Bock et al. 2006), and Spitzer to resolve the "miss- 
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ing source" problem in the near-IR background intensity. We 
conclude with a summary in § 7. Throughout this paper, we 
take a reference cosmology model with parameters Q m = 0.3, 
J7a = 0.7, h = 0.7 (unless left as a free parameter, such as in the 
luminosity function), and a normalization to the matter power 
spectrum at 8 h _1 Mpc scales a% = 0.84. 

2. IMAGING DATA 

To generate the angular power spectrum, we make use 
of the source catalogs in the GOODS and the NDWFS- 
Bootes fields that were extracted using the SExtractor pro- 
gram (Bertin & Arnouts 1996) using 3.6" and 5" apertures, 
respectively. Vega magnitudes were calculated from the cat- 
alog flux values by converting to AB magnitudes and adding 
the appropriate offset from Kashlinsky (2005), and through- 
out this paper magnitudes refer to Vega magnitudes. We refer 
the reader to Table 1 for basic statistics related to our fields 
such as the number of galaxies in specific magnitude ranges 
and the contribution to IRB from those galaxies when the 
counts are integrated. We now provide a brief summary of 
each of the two fields. 

2.1. Bootes field 

We make use of both Spitzer IRAC imaging data as well as 
the ground-based J- and /^-band imaging data of the Bootes 
field from the FLAMEX survey (Elston et al. 2006) where the 
source catalogs are publicly available We refer the reader 
to Eisenhardt et al. (2004) for details related to this IRAC 
Shallow survey and the first results from these data, while 
Fazio et al. (2004a) contains a discussion of number counts. 
The Spitzer Bootes field covers a total area close to 8.5 deg. 2 , 
but due to non-uniform coverage and to avoid introducing a 
complicated window function for clustering measurements in 
the multipole space, we restrict the discussion to the central 
6.3 deg. 2 here. While ground-based FLAMEX imaging data 
of the Bootes field cover a total area of almost 7 deg. 2 , again 
due to gaps in the data, we concentrate on 3.9 deg. 2 with uni- 
form coverage at the center of the field. 

To understand the extent to which Bootes field imaging is 
complete, we performed standard simulations by placing ten 
thousand artificial point sources randomly across the Bootes 
mosaic in each 0.2 magnitude bin. We perform object detec- 
tion and extraction exactly as we had carried out with the orig- 
inal survey catalog as described in Eisenhardt et al. (2004). 
We quantify completeness as the recovered fraction of point 
sources in each magnitude bin. The formal 50% complete- 
ness limit of the catalog occurs at a 3.6^m magnitude of 19.3. 
To the 3.6/im limit of 17.8 mag, the IRAC Shallow survey 
is over 85% complete. Due to completeness and to avoid sys- 
tematic biases from stellar sources in the galaxy that dominate 
at magnitudes below 15 (see, Fazio et al. 2004a) our cluster- 
ing measurements in the Spitzer Bootes field is restricted to 
galaxies in the magnitude range between 17 and 19. In the 
case of ground-based surveys, the 50% completeness limit for 
5 <7 detection of point sources is in K$ band is 19.5 (Elston et 
al. 2006). Thus, we restrict clustering from ground-based im- 
ages to the magnitude range between 17 and 19 in Ks band. 
Since typical galaxy color J—K s is around 1.5 (Elston et al. 
2006), we can extend clustering studies in J-band to magni- 
tudes around 21. In the case of Spitzer L-band galaxy sample, 
typical color L-K is about 1.1 (see, also Oliver et al. 2004). 

10 http://flamingos.astro.ufl.edu/extragalactic/overview.html 



In Figure 1 right panel, we show the completeness corrected 
counts in L-band and those published by Fazio et al. (2004a). 
We find a good agreement with completeness corrected counts 
and those published by Fazio et al. (2004a) using the same 
imaging data. We do not show counts in J and Ks bands as 
these are discussed in detail in Elston et al. (2006). 

2.2. GOODS fields 

Imaging data of GOODS fields (both CDF-S and HDF-N) 
with Spitzer IRAC cover roughly 0.046 deg. 2 . In both fields, 
catalogs of sources in the GOO DS IRAC 3.6/xm mosaics were 
generated using the SExtractor (Be rtin & Arnouts 19961) soft- 
ware. A 9 x 9 Mexican-hat kernel was used with a FWHM of 
2.4" (Dickinson et al. 2006 for further details). Monte-Carlo 
simulations were performed on the IRAC images wherein ar- 
tificial sources, convolved with the Spitzer point spread func- 
tion, were input at random positions on to the mosaic. At 
faint flux densities, almost all sources are point sources at the 
IRAC spatial resolution of 1.6" FWHM (Fazio et al. 2004b). 
100 sources were added in each iteration, which is negligi- 
bly small compared to the ~ 14,000 sources which are in the 
IRAC catalog. The flux distribution of artificial sources was 
flat in log(S3.6) enabling equal numbers of sources in each log- 
arithmic flux bin. The 50completeness corrections outlined in 
the next paragraph, we extend clustering measurements down 
to a magnitude limit of 22.5. 

Sources in the original GOODS catalogs were first matched 
to sources in the catalog generated on the image with the fake 
sources, with a 1" matching threshold. If the sources in the 
original catalog were not matched it implies that the flux and 
position were both affected by confusion/blending from an 
artificial source. The positions of the input artificial sources 
are known and these are matched with the unmatched sources 
in the catalog with the same positional tolerance. An artifi- 
cial source of a particular flux would then be detected if it is 
within 1" of its original position although its extracted flux 
might be significantly different from the input flux if it landed 
close to an existing bright IRAC source. The artificial source 
would, however, be undetected if its position in the catalog is 
more than an arcsecond from the original position. The flux 
distribution of input and output sources is transformed to a 
matrix P,- ; where index i is the input flux and index j is the 
output flux (Chary et al. 2004). The nature of the Py matrix 
is such that for a particular i, the sum over all j is less than 
unity. This is the completeness correction factor that we use 
for the GOODS CDF-S and HDF-N catalogs. 

3. STATISTICAL PROPERTIES OF GOODS AND BOOTES FIELDS 

Here, we make use of two types of surveys: the Bootes field 
provides us with a shallow, wide-field image to measure clus- 
tering out to degree angular scales for bright sources (down 
to a magnitude limit of 19 in the L-band), while the GOODS 
HDF-N and the GOODS CDF-S Spitzer images allow us to 
measure clustering of resolved sources at 10 arcminute angu- 
lar scales and below down to a very deep magnitude limit. 
In combination, the Bootes image allows us to measure clus- 
tering from degree scales to sub-arcminute scales for bright 
galaxies and GOODS catalogs allow clustering measurements 
of resolved sources down to the faintest level of 22.5 in L- 
band. This is at the magnitude level at which resolved sources 
were removed from the fluctuation analysis in Kashlinsky et 
al. (2005). Thus, clustering studies with the two fields, in re- 
turn, allow us to establish the confusion level when making 
measurements related to the unresolved diffuse component 
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FIG. 1 . — Left: Number counts in the GOODS HDF-N and Bootes fields with Spitzer IRAC at L-band, and Bootes fields from the ground at J- and Xj-bands 
using data from the FLAMEX survey, as described in § 2. Right: Comparison to Fazio et al. (2004a) counts and the completeness corrected counts of the Bootes 
field used in the analysis of this paper. Between the magnitude range of 17 and 19 used for clustering measurements in the Bootes field, we find a good agreement 
with counts extracted for the present analysis and the ones studied by Fazio et al. (2004a) for total background intensity measurement at 3.6 fim. Note that Fazio 
et al. (2004a) counts shown here are the completeness corrected ones at magnitudes fainter than 16. 
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FIG. 2. — The fractional (solid symbols; in nW m 2 sr 1 magnitude ') and 
total (open symbols; in nW nT 2 sr~') contribution to the intensity of the IRB 
at 3.6 fim from galaxy counts in GOODS and Bootes fields. The two hori- 
zontal lines show the measured range of the total intensity at 3.6 fim of 12.4 
± 3.2 nW nT 2 sr~' (Wright & Reese 2000). The counts are completeness 
corrected. Using GOODS counts, we measure a total intensity for the IRB at 
3.6 fim of ~ 6.6 nW nT 2 sr~' consistent with previous estimates. 

and to understand if the fluctuations detected by Kashlinsky et 
al. (2005) can be explained as simply extrapolating the faint 
number counts to below the point-source detection threshold 
applied in Kashlinsky et al. (2005). 

3.1. Number counts and the IRB intensity 

In Figure 1 left and right panels, we show the original num- 
ber counts and completeness corrected number counts from 
various catalogs, respectively. In the right panel, the com- 
pleteness corrected counts are compared with those published 
by Fazio et al. (2004a). As shown in Fig. 1 (right panel) 
and for completeness corrected counts in both Bootes and 
GOODS fields, we find good agreement to the extent that we 
can compare the fields reliably. In Fig. 2 we show the frac- 
tional and integrated contribution to the total intensity of the 



IRB in each of the magnitude bins (e.g., Totani et al. 2001). 
For comparison, we also show the total intensity of the IRB at 
3.6 /xm from Wright & Reese (2000). 

Integrating the source counts corrected for the stellar num- 
ber counts, J SdN/dSdS, extragalactic sources between the 
magnitude range of 8 and 17 produce a cumulative back- 
ground of 2.4 nW m" 2 sr" 1 , while incomplete source counts 
between L-band magnitudes of 17 and 19 generate a IRB in- 
tensity of about 1.8 nW m" 2 sr" 1 . The total intensity of Bootes 
source counts between 17 and 19 magnitudes increases to 2.1 
nW m" 2 sr" 1 when counts are corrected for completeness. In 
combination, this is roughly 37% of the total background at 
3.6 /im of about 12 nW m" 2 sr" 1 (Wright & Reese 2000). 

For sources fainter than 19th magnitude, the background 
varies depending on the completeness corrections and the as- 
sumed slope below the point-source detectability level with 
the GOODS images. For completeness corrected source 
counts between L-band magnitudes of 19 and 22.5, we find a 
total background of 2.3 nW m" 2 sr" 1 . If the faint-end GOODS 
slope in Fig. 1 left panel is extrapolated to a magnitude limit 
well below 30, we find a cumulative background of 1 .2 nW 
m" 2 sr" 1 for sources fainter than a magnitude of 22.5. As we 
discuss later in Section 4.2, such a slope is inconsistent with 
fluctuation measurements of the unresolved background by 
Kashlinsky et al. (2005), though such a slope is suggested in 
the study of the one-point distribution function of the unre- 
solved flux by Savage & Oliver (2005). The faint-end slope 
that describes the anisotropy measurements of Kashlinsky et 
al. (2005) suggests a background intensity from sources be- 
low a magnitude limit of 22.5 of 0.08 nW m" 2 sr" 1 . Com- 
pared to the total intensity of about 12 nW m" 2 sr" 1 , we find 
that resolved source counts and unresolved background clus- 
tering explain about 55% of the background for extragalac- 
tic sources, which is consistent with the resolved fraction in 
Fazio et al. (2004a), but lower than the resolved fraction sug- 
gested in Savage & Oliver (2005). 

3.2. Redshift distribution 

The photometric redshift distribution for galaxies with 
magnitudes between 17 and 19 in 3.6/xm is computed as 
described in Brodwin et al. (2006). We show this red- 
shift distribution in Figure 3 left panel with a comparison to 



5 




FIG. 3 . — Left: The normalized (J dzw(z) = 1) redshift distribution of L-band galaxies with magnitudes between 17 and 19 as estimated photometrically for the 
GOODS HDF-N and CDF-S fields (Mobasher et al. 2004) as well as for the Bootes field (Brodwin et al. 2006). We take into account these distributions when 
model fitting the L-band clustering measurement in this magnitude bin. Right: The redshift distribution of ./-band galaxies with magnitudes between 17 and 19 
and between 19 and 21, as estimated again based on photometric redshifts. 



GOODS photometric redshift distributions in the same mag- 
nitude range for both HDF-N and CDF-S fields. The GOODS 
distributions are described in Mobasher et al. (2004) and are 
accurate to (z p hot-Zs P ec)/(l +z sp ec) ~ 0.05 statistically. There 
is a broad agreement between Bootes and GOODS distribu- 
tions with the distributions peaking at redshifts between 0.8 
and 0.9. Based on numerical calculations, we determined that 
the fractional contribution to the clustering power spectrum 
from galaxies at redshifts greater than 2 with magnitudes be- 
tween 17 and 19 in L-band is below a few percent at the angu- 
lar scales probed by Bootes and GOODS catalogs. Thus, with 
a clear peak of the redshift distribution at z ~ 0.9, we can eas- 
ily interpret parameters extracted from galaxy clustering as 
those related to sources at this redshift instead of accounting 
for the corrections coming from the high-redshift tail of the 
redshift distribution. 

Furthermore, while not exact in terms of selection criteria or 
observing wavelengths, the redshift distribution for the bright 
IR galaxy sample between 17 and 19 magnitudes in the L- 
band is similar to the galaxy redshift distribution in the high- 
redshift end of the DEEP2 survey. Thus, our statistical results 
on quantities such as the halo mass hosting IR galaxies can be 
compared with those established based on the DEEP2 clus- 
tering measurements (Coil et al. 2004). In the right panel of 
Figure 3 we show the redshift distributions of /-band galaxies 
with magnitudes between 17 and 19 and between 19 and 21, 
as estimated photometrically. These again peak at z ~ 1 and 
we use them when model fitting the clustering measurements 
in the /-band with source catalogs generated by the FLAMEX 
survey. 

4. CLUSTERING MEASUREMENT TECHNIQUE 

The measurement of galaxy clustering in the multipole 
space is now well described in the literature with a wide range 
of applications with data from several large-scale structure 
surveys. These techniques borrow from ideas that have been 
developed for clustering measurements in the multipole space 
in cosmic microwave background (CMB) anisotropy maps. 
In general, the available techniques fall into two approaches. 
First, starting from a technique described in Peebles (1973), 
clustering can be derived either through spherical harmonic 



estimation with moments calculated by summing over all pix- 
els. The second approach is to apply a maximum likelihood 
method to data to estimate the power spectrum directly, espe- 
cially when data involve large numbers of pixels where direct 
computation of spherical transforms is computationally time 
consuming (e.g., Tegmark et al. 2002). Here, we make use of 
the former approach given the small area of our surveys and 
the ability to directly estimate the transform. 

To begin, in each catalog, the source count distributions 
were fit to a fine grid of pixel size 1.8" for Bootes and 0.7" 
for GOODS, with each pixel i in the grid containing the sum 
of the flux of all the sources contained within it, 7,(n), where 
n is the two-dimensional angular vector on the sky. In com- 
parison to usual clustering measurements of the galaxy distri- 
bution, 7,(n) is not just the number of galaxies in each pixel, 
but rather the total IRB intensity in that pixel given the magni- 
tude of that galaxy. Note that we put all flux within the pixel 
at the source center and do not account for the large IRAC 
PSF While this procedure ignores the spatial structure of indi- 
vidual sources, since we are measuring large-scale clustering 
between sources rather than small-scale clustering at spatial 
scales below a typical source size, this is not a concern for us 
for this study. It would, however, be a concern when studying 
clustering in the unresolved background (Section 6.2). 

In summary, to extract C/, we expand the 7,(n) as 

'max m=l 

;=o m=-l 

where F; m are the spherical harmonic functions. In principle, 
if data exists over the whole sphere, angular power spectrum 
is simply 

i ' 

c ' = ^TYE[l / "»l 2 - c ' hot " noi 1 > & 

m=—l 

where the factor (21 + 1) comes from the fact that at each 
multipole /, there are (21 + 1) independent modes. In above, 
^■shot-noise ^iggg f rom the discrete behavior of the two di- 
mensional intensity field (see below) and acts as a source 
of noise for clustering measurements with a finite density of 
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sources (Peebles 1973) and this shot-noise can be written as 
qhot-noise = j s 2 dN /dSdS, where dN/dS is the surface den- 
sity of sources at flux S, assuming all counts are extracted with 
no confusion. 

In practice, however, one must correct clustering measure- 
ments for the incomplete sky coverage through an appropriate 
renormalization of the harmonics. Following Peebles (1973), 
for partial sky coverage, we can write an estimate for the an- 
gular power spectrum as 



C, 



21 



—Y 
+i ^ 



\hm-C\ 



shot-noise / 



Mm | 



B, 



_ ^shot-noise 



(3) 



m=-l 

where the two functions A/,„ and B\ m are A[ m = J Y[ m (n)dn 
and B[ m = j \Y[ m (h)\ 2 dh, respectively, and the angular inte- 
grals are over the area of the survey. While in the all-sky 
case, estimates of C/ at each multipole / is independent of 
each other, with partial sky coverage, adjacent estimates are 
correlated. Thus, one cannot represent the power spectrum 
with estimates at each I. To decrease the correlations, we 
average adjacent Q estimates over bins Al in the multipole 
space such that Al is taken to be logarithmically dependent 
on central multipole or A = xl, with x is fixed at the level of 
0.2. With such large bins, using the window function in mul- 
tipole space, Wi = Y^m=-i \Aim\ 2 / {11 + 1), we have established 
that the power spectrum estimates are correlated at a level be- 
low a few percent and we neglect these correlations in our 
model fitting analysis. While not pursued here due to simple 
survey geometry and the small area used for clustering mea- 
surements, the above procedure can be optimized further to 
quickly calculate harmonic transforms over large areas on the 
sky using specially designed pixelization schemes that have 
been developed for CMB anisotropy studies (e.g., HEALPix; 
Gorski et al. 2005) and pseudo-C/ estimates (e.g., Wandelt et 
al. 2001). 

As described earlier, to avoid complicating the clustering 
measurement with images that have large gaps in the source 
distribution and a complex window function W/, we chose a 
smaller, central area within the Bootes field of 6.3 square de- 
grees with IRAC catalogs and a 3.9 square degrees field in 
ground-based images. Such a selection also guarantees that 
our clustering measurements are not affected by problems at 
the edge of the field and that we do not need to decorrelate the 
measurements for the complex geometry. For the pixel sizes 
used, based on a Monte Carlo simulation, we have established 
that the contribution to the power spectrum from overlapping 
sources was negligible. However, our grid has empty pixels 
due to masking of bright sources or no data (e.g., muxbleed 
contamination) which we describe through a mask and such 
empty pixels cannot be avoided since they occur throughout 
the field. To avoid introducing complications to the power 
spectrum measurement, we filled those random empty pix- 
els with white noise before estimating the power spectrum, 
where the white noise was estimated using the same proce- 
dure as the boot-strap approach outlined below to establish 
errors for binned C/. To test the extent to which such a pro- 
cedure can affect our clustering measurements, we utilized a 
Monte-Carlo approach by systematically masking out differ- 
ent regions that were not contaminated or, instead of adding 
simple white noise distributed as a Gaussian, replacing the 
contaminated regions with Poisson noise. We found that clus- 
tering measurements do not vary beyond the error indicated 
and that large-scale clustering is unaffected by contaminated 
pixels. 



As outlined above our clustering measurements lead to an 
estimation of the total power spectrum, which includes both 
the shot-noise component, associated with the Poisson-noise 
or discreetness of the source counts, and the clustering com- 
ponent, associated with the true spatial distribution. The lat- 
ter is the component that is of interest to large-scale structure 
studies and, as outlined above, one must remove the shot- 
noise component. There are two approaches to estimate the 
shot-noise. First, one can simply measure the total power 
spectrum including the shot-noise and then account for the 
shot-noise part of the spectrum, c ; sho, ~ nolsf % by measuring the 
amplitude of clustering, that will scale simply as £ 2 when plot- 
ted in terms of / 2 C; /2tt or as a constant is simply considered 
as Q, at large £ values. This assumes that at large £ values or 
small angular scales, clustering is dominated by pure noise. 
Since noise for our measurement is simply the shot-noise, 
this amplitude can be taken to be the shot-noise and that con- 
stant value can be removed from the total spectrum C/. The 
second approach is to use the information related to existing 
number counts and simply do the integral J S 2 dN/dSdS to 
estimate the expected shot-noise, since this is also a correct 
description of the shot-noise. In practice, we have found that 
the shot-noise component is more accurately measured using 
the amplitude of the total power spectrum at smallest angular 
scales probed by the data just before the finite size of pixels 
become important than simply evaluating the integral over the 
number counts in the magnitude range for which clustering 
measurements are made. This is due to the incompleteness in 
our measured number counts as we now describe. 

To understand this, our clustering spectrum is estimated 
from the IRB intensity captured by individual sources with 
the image pixelized in terms of flux. Incompleteness in the 
point source detection catalog leads to a systematic under- 
estimate of the flux in certain pixels. In addition to incom- 
pleteness, our clustering spectrum is also affected by stellar 
sources that were not removed in the original catalog. In this 
case we may be overestimating the clustering amplitude from 
contamination related to stellar flux. Thus, in addition to in- 
completeness of the sky coverage which is easy to correct for 
once the geometry of the survey area is specified, we use a 
two-step approach to correct for the angular power spectrum 
of galaxies in each of the magnitude bins. First, our clustering 
measurements are corrected to account for incompleteness in 
the catalogs. Using the difference between J S 2 dN/dSdS for 
completeness corrected counts and the one estimated directly 
from clustering data that is an underestimate of the total shot- 
noise, we correct the shot-noise and clustering power-spectra 
as 

/^shot-noise _ /.shot-noise / S 2 dN / dS COmplete dS 

-<-/ 



Q = Ci 



J S 2 dN/dS incom v lele dS 
J SdN/dS com v lete dS 



J SdN/dS incom P lete dS 



(4) 



When making these corrections to the overall amplitude we 
are assuming that incompleteness affects the field uniformly 
and that incompleteness correction has no angular structure. 
In optical surveys, incompleteness can be related to extinction 
and for large are surveys such as SDSS that span over 1000 
deg. 2 or more, extinction varies across the field significantly 
and these angular variations must be accounted for. For the 
present study involving IR galaxy clustering at angular scales 
of a few degree and below, such complex corrections can be 
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safely ignored. 

The completeness corrected counts in Figure 1 left panel 
results in a 18% increase in the shot-noise, which if not ac- 
counted for leads to an underestimate of power at small an- 
gular scales. At arcminute scales where clustering measure- 
ments from the Bootes field overlap with the GOODS, we find 
good agreement between the measurements from these two 
fields when Bootes clustering measurements are corrected for 
completeness issues. While it is not necessary to correct for 
completeness in the GOODS fields between magnitudes of 17 
and 19 in the L-band, for clustering measurements in fainter 
bins, we make use of the completeness calculations to correct 
for the shot-noise and the amplitude. 

Since we did not correct the catalogs for potential contam- 
ination from stars, the shot-noise contains an additional con- 
tribution from the stellar counts in addition to galaxies whose 
clustering we are attempting to measure. To correct for the 
shot-noise associated with stars, we make use of the model 
of Fazio et al. (2004a) and calculate the shot-noise associ- 
ated with stars in the magnitude range of 17 to 19. The tabu- 
lated counts in Fazio et al. (2004a), based on two independent 
deep fields, lead to two different estimates of the stellar shot- 
noise correction. The stellar shot-noise is 6% and 10% of 
the original total shot-noise for the two fields and is a small, 
though non-negligible, correction to the angular power spec- 
trum. We show both cases when presenting our clustering 
measurements to highlight the small difference, which is only 
important at arcminute and smaller angular scales. Note that 
the corrections related to stellar sources are done similar to 
that of the catalog incompleteness and again we assume that 
confusion from stars are uniform throughout the image. 

In Figure 4, we show only the clustering component as this 
is directly related to the large-scale structure of the source 
distribution. The measurements are binned in multipole. Er- 
rors in each bin were computed through a boot-strap method, 
where a part of the source distribution was replaced by white 
noise and the power spectrum measured 100 times by repeat- 
ing the process. The power spectrum error is estimated by the 
bootstrap sample variance associated with the multiple mea- 
surements and we have established that these errors are con- 
sistent with the theoretically expected error for C; of 

to better than 10% when / sky is the fraction of sky covered in 
the survey and Al is the bin width for each C; measurement. 

To test how well we have estimated errors in our clustering 
spectrum, in Figure 5 we plot the rms fluctuation in the IRB 
from resolved counts with magnitudes between 17 and 19 at 
3.6^m at an angular scale of ~ 4.5 arcminutes (£ = 7200; left 
panel) and ~ 1.2 arcminutes (I = 28800; right panel) when 
measured in 100 independent fields of 0.05 degrees randomly 
extracted from the Spitzer IRAC Bootes image. For compar- 
ison, we also show the rms fluctuation level estimate at these 
multipoles from GOODS HDF-N and the estimated la error 
bar of the rms based on the jackknife sampling described in 
the text. The variance in 100 independent fields within Bootes 
of GOODS-sized fields are consistent with the error estimated 
from GOODS data, suggesting that we are properly account- 
ing for the cosmic variance. This comparison suggests that 
our method to estimate errors leads to reasonable estimates of 
the uncertainties in the clustering spectrum. 

In Table 1, we list the magnitude bins considered for clus- 



tering measurements with the GOODS HDF-N and the Bootes 
catalogs at L-band from Spitzer IRAC and J- and /Tj-bands 
from ground-based imaging of the Bootes field. The first 
bin is selected to provide an overlap between the two cata- 
logs. Bootes with the largest sky coverage of the surveys al- 
lows large-scale linear clustering measurements to be made, 
while GOODS allows measurements in the mildly non-linear 
to non-linear regime. In Table 1, we also list the mean IR 
background produced by sources in each of the bins as well 
as the shot-noise to the anisotropy measurements; the for- 
mer comes from the source count distributions while the latter 
comes from direct anisotropy measurements at small angular 
scales. 

5. ANALYTICAL MODEL FOR IR SOURCE CLUSTERING 

Following Cooray et al. (2004), we consider a simple 
halo model-based description for the angular clustering of IR 
sources. The angular power spectrum can be written as 

«-/*s^*w-(»-ir")- ,6) 

where a(z) = (1 +z)~ l is the scale factor, cIa is the comoving an- 
gular diameter distance and r is the radial distance, while i\(z) 
is the mean emissivity per comoving volume of IR sources at 
wavelength A and P ss (k) is the power spectrum of the clus- 
tered sources. Based on a comparison to the exact formula 
involving integrals over spherical Bessel functions, the Lim- 
ber approximation (Limber 1954) used here to describe the 
clustering spectrum of IR sources in the multipole space is 
accurate to better than 0.2% at multipoles less than 10 and is 
accurate to even higher precision for degree scale fluctuations 
at i ~ 10 : . 

Instead of modeling the exact relation between IR lumi- 
nosity and halo mass, we assume that over a given magni- 
tude range, we can simply scale source fluctuations by the 
mean IR background produced by those sources to get back 
to fluctuations in the IR intensity. In this case, the angular 
power spectrum takes a simpler form of C/ = Pwi, where 
P is the cumulative intensity of sources in the magnitude 
range of interest and wi is the projected angular power spec- 
trum of the sources relative to the density field such that 

wi = J drw 2 (z) J 'd\P ss (k= -^,zj, where w(z) is the normal- 
ized radial distribution of sources such that J dzw(z) = 1 . We 
make use of the statistical redshift distribution implied by the 
photometric redshift estimates for galaxies in the Bootes cat- 
alog (see, Figure 3). 

In addition to the halo model, clustering can also be 
predicted based on linear clustering spectrum that ignores 
growth of perturbations and non-linearities that are gen- 
erated at low redshifts as structures grow under gravita- 
tional evolution. For this, we simply take the linear spec- 
trum in our fiducial cosmology, and then calculate C ; lln = 

P jdrw 2 (z)/d 2 P M (k=j I ,z^, where D 2 (z)P M (k) is the lin- 
ear power spectrum and D(z) is the linear theory growth func- 
tion (Peebles 1980). 

This linear theory spectrum, however, is not an appropriate 
description since it fails to describe clustering at small angu- 
lar scales or multipoles above a few hundred. Thus, more 
accurate descriptions of the underlying source spectrum is re- 
quired. For this, we make use of the halo model where the 
three-dimensional source power spectrum contains two terms 
with P ss (k) = P lh (k) + P 2h (k) (see, Cooray & Sheth 2002). 
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FIG. 4. — Left: L-band clustering measurements in both the GOODS fields and the NDWFS Bootes field with Spitzer IRAC. The measurements are divided 
into magnitude bins in the catalog from 17 to 19 (45 to 7 /ijy) (for both Bootes and GOODS; top measurements), 19 to 21 (7 to 1.1 /ijy) (GOODS; middle), and 
21 to 22.5 (1.1 to 0.28 /^Jy) (GOODS; bottom). The angular scale of clustering is related to the multipole I through roughly 8 ~ n/£. In the case of GOODS, 
filled symbols are the clustering measurements in the HDF-N while open symbols are from the CDF-S. The Bootes clustering spectrum is extracted based on 
the completeness corrected shot-noise following counts in Figure 1 left panel. The open and filled symbols associated with the Bootes clustering spectrum show 
the difference between two estimates for the stellar contribution to the number counts and the shot-noise. This contribution is 6% of the total shot-noise in the 
magnitude range of 17 to 19 in L-band for the open symbols and 10% for the filled symbols, and this difference is relatively insignificant. The lines show model 
descriptions of the angular power spectrum of source clustering using CLF models in 17 to 19 and 21 to 22.5 magnitude bins (see Section 5.2) and we show the 
total (solid lines), 1-halo term (dot-dashed lines), and the prediction under linear theory for galaxy clustering (dashed lines). Right: Clustering of resolved point 
sources in the ground-based imaging data of the Bootes field at J and Ks bands in magnitude bins. Note the factor of 100 difference in the y-axis scale between 
the left and the right panels. This difference is consistent with the mean intensity of the IR background between these two different wavelengths. The curves here 
show the total (solid and dot-dashed for J and Ks-band clustering between 17 and 19 magnitudes, respectively), and for reference, in the case of J-band clustering, 
the expected linear model prediction. 
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FIG. 5 . — The rms fluctuation in the IRB from resolved counts with magnitudes between 17 and 19 at 3.6pm at an angular scale of ~ 4.5 arcminutes (£ = 7200; 
left panel) and ~ 1.2 arcminutes (i = 28800; right panel) when measured in 100 independent fields of ~ 0.05 degrees randomly extracted from the Spitzer IRAC 
Bootes image. For comparison, we also show the rms fluctuation level estimate at these multipoles from GOODS HDF-N and the estimated 1 cr error bar of the 
rms based on the boot-strap sampling described in the text. The variance in 100 independent fields within Bootes of GOODS-sized fields are consistent with 
the single error estimated for the GOODS HDF-N catalog directly through bootstrap sampling. This agreement suggests that we are properly accounting for the 
cosmic variance of the angular power spectrum. 



These two terms are clustering between IR galaxies in two 
different halos (2h) and clustering of galaxies within the same 
halo (lh), and given by 



P 2h (k) = 



dM n(M) b(M)^^—^-u(k\M) 



P lin (k) 



,„ / 2(N,)(N c )u(k\M)+(N s ) 2 u 2 (k\M) 
P lh {k)= / dMn(M) V A 1 V 1 ' V ' — — ,(7) 



respectively. Here, u(k\M) is the normalized density profile 
in Fourier space (e.g., Navarro, Frenk & White 1996), n(M) 
is the halo mass function (e.g., PS mass function of Press & 
Schechter 1974), b(M) is the halo bias relative to the linear 
density field (e.g., Mo, Jing & White 1997), and n g is the 
number densities of galaxies defined below. As written, the 
2-halo term with P 2h (k) traces the linear power spectrum, but 
at large k, through u(k\M), it departs from the linear theory 
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TABLE 1 



Survey 


Band 


Area 
[deg. 2 ] 


Mag. Range 
[Vega] 


N 

1 * sources 


A 1(A) 
[nWirf 2 sr _1 ] 


/ S 2 dN/dSdS 
[nW 2 nT 4 sr _1 ] 


[A 1(A)]'' 
[nWnT 2 sr 1 ] 


[ / S 2 dN/dSdS] c 
[nW 2 nT 4 sr 1 ] 


NDWFS 


L 


6.3 


17-19 


2.2 xlO 5 


1.75 


3.36 x 10~ 8 


2.09 


3.96 x 10~ 8 


(Bootes) 


K s 


3.9 


17-19 


6.3 xlO 4 


2.9 


2.05 xlO" 7 








./ 


3.9 


17-19 


2.7 x 10 4 


5.51 


1.72 xl0~ 6 








./ 


3.9 


19-21 


6.2 x 10 4 


2.13 


1.10 xlO~ 7 






GOODS 


L 


0.046 


17-19 


2.1 xlO 3 


1.99 


3.69 x 10~ 8 


2.20 


4.06 x 10~ 8 


(HDF-N) 


L 


0.046 


19-21 


5.5 x 10 3 


0.86 


2.65 x 10~ 9 


1.24 


3.69 x 10"' 




L 


0.046 


21-22.5 


4.8 xlO 3 


0.17 


1.02 xlO~ 10 


0.39 


2.26 x 10~ 10 



NOTES. — Survey parameters associated with clustering measurement in both GOODS (HDFN; Dickinson et al. 2003) 
and NDWFS Bootes (Jannuzi & Dey 1999) fields in J-, K s , and L-bands. With deep GOODS HDF-N data we consider 
three bins in magnitudes in the L-band. The L-band Bootes field data are from Spitzer while /- and Ks data involve 
catalogs of sources detected with ground-based imaging of the Bootes field by the FLAMINGOS Extragalactic Survey 
(FLAMEX; Elston et al. 2006). The last four columns tabulate the mean background intensity and the shot-noise for 
uncorrected and corrected counts (with [A I(A)] r and [intS 2 dN /dSdSf ), respectively (see, Section 3.1 for details on 
completeness corrections). 
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FIG. 6. — 3.6 fim galaxy luminosity functions (LFs) as a function of 
redshift. Plotted data points are from Babbedge et al. (2006). The solid 
curves show the LF based on a conditional luminosity function model. For 
< z < 0.25, we also show the LF from central (dashed line) and satellite 
(dot-dashed line) separately. 

prediction alone. 

To calculate halo occupation numbers, (A^ s ) and (N c ) with 
(N ga i(M)) = (N s ) + (N c ), we consider two separate ways to 
model the galaxy distribution. First, we consider halo occupa- 
tion numbers where all galaxies are treated the same regard- 
less of their luminosity or flux and then built a more com- 
plete model using conditional luminosity functions (Cooray 
2006a) that are normalized to reproduce the 3.6 /im galaxy 
LFs measured in Spitzer SWIRE surveys out to a redshift of 
1.0 (Babbedge etal. 2006). 

5.1. Simple halo model 

In the simplest description for the halo occupation, the 
source distribution within halos is described analytically as 
(N ga i(M)) = 1 + (N S (M)) whenM > M min and otherwise, with 
the assumption of a central galaxy in each halo ((N C (M)) = 1) 



and a power-law distribution of satellites that scales with halo 
mass as (N S (M)) = A(M/M min ) /3 ; We take parameters A, (3, and 
M m i n to be free parameters to be determined from the data. 
Here, M m i n is the minimum halo mass scale at which galax- 
ies appear over the magnitude range of interest; this is not 
the same as the minimum mass scale at which any galaxy can 
appear since smaller halos may host galaxies with lower lu- 
minosities, but we do expect an overall cut-off for the appear- 
ance of any galaxy at some low mass scale for dark matter 
halos (see Section 6.2). The mean number density of galaxies 
is fig = J dM n(M) (N ga i(M)} . 



5.2. Conditional luminosity functions 

While the above description does not separate galaxies 
based on their luminosity, we can extend our formalism to 
separate the occupation number to a fixed central galaxy and 
satellites around them through conditional luminosity func- 
tion models (CLF; Cooray & Milosavljevic 2005; Cooray 
2006a,b; van den Bosch et al. 2006). These CLFs describe 
the average number of galaxies with luminosities between L 
and L + dL that resides in halos of mass m at a redshift of 
z, while halo occupation numbers defined above describe the 
total number of galaxies in a dark matter halo of mass m re- 
gardless of the luminosity of those galaxies. The CLF models 
are designed to model not just clustering statistics, but also 
one-point moments of the distribution such as LFs and num- 
ber counts. We refer the reader to Cooray (2006a) for details, 
but here provide a summary and related fitting functions. 

We define redshift-dependent CLFs such that 



$(L|M,z) = $ c (L|M,z) + $ s (L|M,z) 
$ c (L|M,z) = ' 



exp 

$ s (L|M,z)=A(M,z)L 7<M) 



27iTn(10)<7 rmc L 

log 1Q [L/L c (M,z)] 2 

2(7 c 



(8) 
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Given the CLF, the galaxy LF is obtained through 

$(L,z) = JdM^(z)HL\M,z). (9) 

Finally, given the information related to observed magnitude 
ranges, one can integrate over the corresponding luminosity 
distribution at a given redshift to obtain a more accurate esti- 
mate of the occupation numbers 

(N c (M,z))= [ dL$ c (L\M,z) 



(10) 



(N s (M,z)) = J dL$ s (L\M,z). 



These can then replace the calculation procedure outlined at 
the beginning of this Section to calculate the angular power 
spectrum. In above, central galaxy CLF takes a log-normal 
shape while satellite CLF takes a power-law behavior with lu- 
minosity. Note that the above two integrals over a log-normal 
and a power-law functions in luminosity, generally recover 
the mass dependence of the simple halo occupation number 
written in the previous section with (N g!l \(M)} = 1 + (N S (M)) 
where 1 is for the central galaxy with (N c (M,z)} and (N S (M)) 
is a power-law with mass. 

An important ingredient in CLF models is the mapping be- 
tween the central galaxy luminosity L c and halo mass (Vale & 
Ostriker 2004) as a function of redshift (Cooray 2006b). This 
relation is described with a general fitting formula given by 



L c (M,z) = L (l+z) a 



{MjM x ) a 



[b+(M/M l y ( ~ l+ ^] l / d 



(11) 



Following the procedures outlined in Cooray (2006a,b), we 
first model the LF of 3.6 /im galaxies between redshifts of 
and 1 (from SWIRE survey; Babbedge et al. 2006). In Fig- 
ure 6, we show the measured LF and the model based on CLFs 
with parameters for the L c (M,z) relation of Lq = 1 .3 x 10 9 L Q , 
Mi = 10 n M Q , a = 4.0, b = 0.7, c = 3.75, and d = 0.22. The 
redshift dependence is captured by fitting parameters a and 
(3 (see, Cooray 2006b) and at 3.6 /im we have found these 
to take values of -0.2 and -0.05, respectively, while the dis- 
persion of central galaxy luminosity (at a fixed mass M) in 
log-luminosity units as written in equation 8 is a = 0.2. As 
discussed in Cooray (2006b) there are large degeneracies be- 
tween these parameters when fitting to the LF and, to obtain 
an overall analytical model to compare with clustering mea- 
surements, we only consider the best-fit description with LF 
curves shown in Fig. 6. 

Finally, we also need a description of the total luminosity 
content of a halo, L tot (M, z), and this is done by using the same 
fitting function as for central galaxies in equation [TT] with 
c = 3.5 while all parameters are kept the same. The luminosity 
content of satellites is then L s (M,z) = L tot (M,z)-L c (M,z), and 
this fixes the normalization A(M,z) in equation[8]by requiring 

that L s (M,z) = J^"™* $ S!it (L\M,z)LdL. When written the satel- 
lite CLF, we take 7 to be independent of mass with a value 
of -1. This is consistent with optical cluster LFs, but these 
models can be improved significantly once similar LF mea- 
surements of galaxy groups and clusters are made at Spitzer 
bands and not just the field LF that averages over a large dis- 
tribution of halo masses. 

6. RESULTS & DISCUSSION 



6.1. Clustering 

In Figure 4 (left panel), we show the clustering measure- 
ments and a model description of the data. GOODS data, 
due to low sky area of 0.092 deg. 2 between the two HDF- 
N and CDF-S fields, measure only non-linear clustering, but 
between magnitudes 17 and 19 (45 and 7 /iJy), Bootes data 
allow the angular power spectrum to be extended to large an- 
gular scales related to the linear part of the angular power 
spectrum. The combination, Bootes and GOODS, shows ev- 
idence for the transition between linear to non-linear cluster- 
ing with a departure from a single power-law clustering spec- 
trum. With counts corrected for completeness, over the angu- 
lar scales of overlap, we find a good agreement between the 
clustering spectrum determined in the large-area Bootes field 
and the ones determined with GOODS images. If not com- 
pleteness corrected, we find a ~ 10% offset between Bootes 
clustering spectrum and the GOODS clustering spectrum at 
multipoles of 10 4 , with Bootes field underestimating power 
relative to the GOODS measurements. 

In addition to the combination of GOODS and Bootes in 
the L-band, we also detect a departure from a power-law in 
the ground-based FLAMEX Bootes /-band data though /ri- 
band clustering measurements could be fitted with a power- 
law (Figure 4 right panel). In Figure 4, the lines show model 
descriptions of the angular power spectrum of source cluster- 
ing using the CLF model for 3.6 /im galaxies in the bin of 17 
to 19 magnitudes and 21 to 21.5 magnitudes separately. These 
predictions make use of the statistical redshift distribution for 
L-band galaxies estimated based on photometric redshift esti- 
mates shown in Figure 3 (left panel). Our model fits are not 
strongly sensitive to assumptions about the redshift distribu- 
tion as long as we do not take all sources to be either at very 
low or very high redshift. Note the factor of 100 difference in 
the y-axis scale between the left and the right panels of Fig- 
ure 4; for sources in the same magnitude bin of 17 to 19 (255 
to 40 /iJy /-band, and 107 to 17 /iJy Z^-band), /-band sources 
produce a factor of ~ 7.5 larger rms fluctuations, y / / 2 C//27r, 
in the resolved IRB relative to sources in the same magnitude 
ranges in the L-band. This large increase is consistent with 
the difference in the mean intensity of the IR background be- 
tween these two different wavelengths and differences in the 
average clustering bias of /-band and L-band galaxies. 

The CLF model-based descriptions of IR galaxy clustering 
spectra, as shown in Figure 4, require knowledge on the re- 
lation between both central galaxy luminosity and halo mass 
and the total luminosity and halo mass. As outlined in Sec- 
tion 5.2, we obtain these relations based on model descrip- 
tions to the measured LFs at 3.6 /im in the Spitzer SWIRE 
survey (Babbedge et al. 2006). This procedure is similar to 
the ones used to described galaxy clustering in SDSS at low 
redshifts and DEEP-2 and Lyman-break galaxies at high red- 
shifts (Cooray 2006a; Cooray & Ouchi 2006). A different 
modeling technique, but with the same underlying approach, 
is available in van den Bosch et al (2006) where CLFs are 
a priori assumed to be Schechter-function shapes. Similarly, 
we also use the simple halo model (in Section 5.1) to model 
fit the data as well. Here, the models require three parame- 
ters, M m i n , the minimum dark matter halo mass at which IR 
galaxies begin to appear, (3, the power-law slope of the satel- 
lite occupation with halo mass, and A, the normalization of 
the satellite occupation number relative to the central galaxy 
occupation of unity in halos with mass above M m ; n . Instead 
of model fits to the CLFs, which involve a large number of 
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parameters, here we will directly constrain parameters such 
as A and (3. The right panel of Figure 4 shows clustering mea- 
surements at /- and /fj-bands using the source catalogs from 
the ground-based imaging data of the Bootes field. Again, 
in addition to clustering measurements, we also show model 
descriptions based on the halo model for galaxies with mag- 
nitude brighter than 19 in both /- and Xj-bands. We have not 
attempted to model fit fainter galaxy samples due to issues 
related to the completeness of the catalogs. 

In Figure 7, we consider a likelihood model fit to the angu- 
lar clustering measurements only by taking two parameters, (3 
involving the slope of the satellite counts as a function of halo 
mass and A the overall normalization of the central-to-satellite 
galaxy occupation number to be free while M m \ n is taken to 
be fixed at 5 x 10 11 M Q for galaxies with magnitudes in the 
range 17 to 19. This M m i n value is not the minimum mass 
for galaxies with any luminosity to appear in halos but rather 
those in the appropriate magnitude ranges. Since we are con- 
sidering IR bright galaxies in the magnitude range between 
17 and 19, this minimum mass is consistent with the expected 
value with the CLF model, though the cut off is not sharp at 
a fixed mass scale given the dispersion we allow between lu- 
minosity and halo mass . With large catalogs, the accuracy of 
clustering measurements can be improved and further param- 
eters determined, but at this stage, it is unlikely that more than 
2 parameters can be extracted from the data. We only model 
fit the range of 17 to 19 magnitude in both the L-band and the 
/-band, since at deeper bins with either the GOODS data or 
ground-based Bootes data alone, the clustering measurements 
are mostly limited to non-linear scales only leading to large 
degeneracies between parameters in the occupation number, 
and the photometric estimates for the source distribution of 
these galaxies are less certain than bright galaxies in the mag- 
nitude bin between 17 and 19. 

Based on the constraints on A and (3 shown in Figure 7, we 
estimate the halo mass scale at which two IR sources with 
a magnitude between 17 and 19 in the L-band appear in the 
same halo, or when (N gal (m)} = 2, is (0.9-7) x I0 li M Q . This 
is consistent with the number in the CLF-based model with a 
value of <~ 2 x 10 12 M Q . Note that in CLF-based models, due 
to the large number of parameters involved, we do not vary 
them to fit them individually, but obtain an overall descrip- 
tion consistent with 3.6 /im LFs between redshifts of and 
1 and also use a single set of values for the input parameters 
when comparing to measurements. In terms of CLF models, 
this is the best we can do at this stage. The mass scale we 
estimate for galaxies at z of 1 can be compared to previous 
estimates in the literature. For example, using conditional lu- 
minosity function based models of galaxy clustering (Cooray 
2006a) compared to DEEP2 (Coil et al. 2004) and COMBO- 
17 (Phleps et al. 2003) measurements of galaxy clustering 
at redshifts between 0.7 and 1.3 suggests that the halo mass 
scale at which bright galaxies with luminosities above L* at 
Z ~ 1 start to appear as satellites is <~ (1 — 8) x 10 12 M Q at the 
lcr confidence range, in good agreement with the same scale 
we have established forz <~ 1 IR galaxies. 

The satellite halo occupation number for IR sources in the 
L-band is such that the power-law slope with mass is below 
0.98 at the lcr confidence level. For /-band galaxies with 
magnitudes between 17 and 19, the power-law slope with 
mass of the satellite occupation number is 0.27 < (3 < 0.89 
at the la confidence level, which is consistent with the ex- 
pected value for the slope 0.72 with halo mass from out 
CLFs. As shown in Figure 7 right panel, we illustrate the 



degeneracy direction in the (3-A parameter plane by plotting 
lines of constant average density of galaxies calculated as 
n g = J dzw(z) J dMdN/dM(N g - d \), where w(z) is the normal- 
ized redshift distribution of sources used for clustering mea- 
surements. The degeneracy in (3-A is such that one traces con- 
stant n g values. This number density, as a function of redshift, 
is directly available with LFs and was the basis to normal- 
ize our CLF description. In future, if LF measurements can 
be improved (the measurements have low signal-to-noise at 
z > 0.5 as can be seen in Fig. 5) and also be extended to higher 
redshifts, then we can hope to improve constraints on the oc- 
cupation number, or more importantly on the conditional lu- 
minosity function. An analysis based on such an approach for 
z ~ 3 to 4 Lyman-break galaxies are available in Cooray & 
Ouchi (2006) and we hope to return to this topic again in a 
later paper using additional measurements we have now be- 
gun to make from the same set of data as used for the cluster- 
ing measurements here (Cooray et al., in preparation). 

6.2. Anisotropics in the IRB 

To understand how these fluctuations of resolved sources 
compare with that of the diffuse IRB, in Figure 8 left panel 
we compare our clustering measurements with the anisotropy 
power spectrum presented in Kashlinsky et al. (2005) for 
the unresolved component with all sources brighter than 22.5 
magnitude (0.28 ^Jy) removed. In the case of these unre- 
solved anisotropics, the small scale structure is that of a shot- 
noise with the power spectrum scaling as £ 2 ; In fact, it is this 
shot-noise that Kashlinsky et al. (2005) used to determine 
that they had removed resolved sources down to a magnitude 
limit of 22.5. The flattening of the fluctuations at the two low- 
est I bins is taken to be indicative of an excess component in 
the IRB, with that component ascribed to first galaxies. The 
predictions shown in lines come from the CLF model with 
parameters that reproduce the LF. 

The right panel of Figure 8 shows the completeness cor- 
rected number counts. The lines here are expectation for num- 
ber counts using the CLF by integrating over LFs as a func- 
tion of volume element. Note that at magnitudes below 18, 
the models slightly under predict the counts despite a reason- 
able agreement with LFs. This difference, however, comes 
from stellar contamination to source counts as our counts are 
not corrected for stars, though based on the stellar model of 
Fazio et al. (2004), we do correct for stars in our clustering 
measurements. 

We show three specific CLF models for counts which are 
all normalized to LFs at z < 1, but with differences at the 
faint end not probed by the LFs. The corresponding clus- 
tering spectra associated with those counts are shown in the 
left panel of Fig. 8. In Fig. 8 right panel, the black lines in- 
volve a model where we set the minimum mass for appear- 
ance of at least one galaxy regardless of the luminosity in a 
dark matter halo to be 3 x 10 10 M Q (or effectively a lumi- 
nosity with a corresponding magnitude below 30 in L-band 
since our discussion below is for galaxies brighter than 30th 
magnitude). That is, at halo masses below this value galax- 
ies are not expected to form due to a variety of reasons from 
photoionization during reionization, preheating of gas from 
AGNs, and/or expulsion of gas from first generation of su- 
pernovae. This number is consistent with the value generally 
quoted in the literature based on semi-analytical models of 
galaxy formation (e.g., Kauffmann et al. 1993; Benson et 
al. 2001). Furthermore, this value is constrained to be below 
3 x 10 10 M Q using direct likelihood fits to clustering measure- 
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FIG. 7. — (a) Model fits to the L-band Bootes and GOODS HDF-N clustering measurements in the magnitude bin of 17 to 19 (45 to 7.0 fiJy), and (b) /-band 
Bootes field clustering in the magnitude bin of 17 to 19 (255 to 40 fi3y). We show the constraints on parameters that describe the satellite occupation number, 
following Section 5.1 (see text for details). In (a), for the L-band galaxies with magnitudes between 17 and 19, model fits suggest fi < 0.98 at the \a level, while 
A is between 0.4 and 1.5. In (b), for /-band galaxies with magnitudes between 17 and 19, 0.27 < f3 < 0.89 at the la level, while A is between 0.4 and 1.1. 
In (b), the dashed lines show constant n s = J dz.w(z) f dN /dM{N sd \) / J dzw(z)\ the degeneracy direction in j3-A plane traces the number density of IR sources, 
normalized to the redshift distribution of the sources, w(z). 



ments in SDSS (Cooray 2006a), though this is for galaxies in 
the r'-band. With this cut-off in halo mass, number counts 
peak at L-band magnitudes around 24. 

To obtain an excess of galaxies at the faint-end, we arbi- 
trarily lower this low mass cut off for galaxy formation and 
set it at a mass scale of 10 8 M Q and take the central galaxy 
luminosity-halo mass relation to continue down to such low 
masses. This introduces an excess of faint galaxies with an 
increase in the number counts at magnitudes starting from a 
magnitude of 22. If we use these counts down to a magni- 
tude limit of 26, the sources between the magnitude range of 
22.5 and 26 produce a shot-noise of ~ 5 x 10" 11 nW 2 m" 4 
sr" 1 , which is roughly a factor of 8 higher than the shot-noise 
seen in the fluctuation measurements (~ 7.2 x 10~ 12 nW 2 m" 4 
sr" 1 ). Note that such a large slope at the faint-end has been 
suggested based on a study of the one-point probability dis- 
tribution function of the unresolved IR intensity (Savage & 
Oliver 2005), the large difference between the expected and 
measured shot-noise level suggests that L-band source counts 
at the faint-end do not continue to increase down to magni- 
tude limits of 26 with the same slope as the one suggested 
with completeness-corrected counts between 19 and 22. The 
order of magnitude lower shot-noise measured by Kashlinsky 
et al. (2005) suggests that the counts either flatten or turn over 
consistent with our basic CLF model with a minimum mass 
cut off around 3 x 10 10 M Q for appearance of IR galaxies in 
dark matter halos. 

The integrated background associated with the default de- 
scription for number counts (in black lines) is ~ 0.08 nW m" 2 
sr" 1 between the magnitude range of 22.5 and 30, suggest- 
ing that if the excess clustering suggested by Kashlinsky et 
al. (2005) is correct, then these sources are unlikely to be a 
significant fraction of the total IRB intensity at 3.6 microns. 
If sources were to continue with the large positive slope to the 
faint-end, then the background between magnitude limits of 
22.5 and 30 would have been close to 1.2 nW m" 2 sr" 1 , which 
is an appreciable fraction given that the total background is 
about 10 nW m" 2 sr" 1 at 3.6 /im. Note that the background in- 
tensity of ~ 0.08 nW m" 2 sr" 1 we suggest for faint sources to 



describe unresolved source clustering shot-noise is also below 
the level of 1 to 2 nW m" 2 sr" 1 background intensity from a 
population of z > 5 galaxies that was considered to reproduce 
clustering measurements of Kashlinsky et al. (2005) by Sal- 
vaterra et al. (2006). Salvaterra et al. (2006) models assume 
that shot-noise is dominating from £ of 10 5 and only consider 
the contribution coming from galaxies at high redshifts. 

To emphasize the redshift ranges where the contributions 
are coming from, in our CLF, we put a low redshift cut at 
z ~ 5. The expected counts are show in green in the right 
panel with corresponding clustering spectrum for those galax- 
ies in the left panel of Figure 8. The clustering level is well 
below the measured level at two orders of magnitude below 
the amplitude of Kashlinsky et al. (2005). Since the cluster- 
ing spectrum is proportional to square of the background, see 
Section 5, one can arbitrarily renormalize the counts, as in 
the case of models in Salvaterra et al. (2006), to fit Kashlin- 
sky et al. (2005) data, but then the observed excess lead to 
a total intensity of IRB an order of magnitude higher. Simi- 
larly, counts are also increased by an order of magnitude. This 
renormaization, however, is not correct since it leads to a large 
surface density of galaxies at redshifts more than 5 at a level 
more than allowed by the expected density of dark matter ha- 
los at these redshifts. This expected high density of galax- 
ies has been already realized in Salvaterra & Ferrara (2006) 
who have argued against a large density of first galaxies that 
dominate the IR background given constraints on the number 
counts. This arbitrary normalization is not an issues for our 
CLF models since the galaxy density is constrained given the 
connection between dark matter halos and galaxies they host. 
Unlike the simple model in Salvaterra et al. (2006) the ampli- 
tude is not a free parameter. 

More importantly, Figure 8 left panel also shows that the 
clustering predicted for galaxies in the magnitude range be- 
tween 22.5 and 26 is consistent with measurements of Kash- 
linsky et al. (2005) for unresolved background and for cluster- 
ing of resolved galaxies. Here, non-linear clustering of galax- 
ies within dark matter halos continue to be important out to 
multipoles of 10 6 . For comparison, in Figure 8 left panel, we 
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also show the clustering measurements in our faintest bin be- 
tween 21 and 22.5 in GOODS. The excess seen in Kashlinsky 
et al. (2005) can be described as simply an extension of re- 
solved component that we have measured and extending to a 
larger angular scale than for measurements in Kashlinsky et 
al. (2005). In combination, there is a consistent description to 
the data from our CLFs suggesting that the unresolved com- 
ponent seen by Kashlinsky et al. (2005) can be described by 
galaxy counts with magnitudes between 22.5 and 26 at red- 
shifts below 5 with contributions mostly from redshifts 1 to 3, 
and that there is no need to invoke a large surface density of 
galaxies at redshifts above 5. 

The limited range of multipoles probed by Kashlinsky et al. 
(2005) limits us from making an exact statement given that 
the data can also be modeled from a clustering spectrum that 
trace the linear density field at redshifts above 8 (dot-dashed 
lines following Cooray et al. 2004). This diffuse background 
is taken to be produced by Pop-Ill stars and their spectra are 
completely different from known galaxies given the dominant 
energy in the Lyman-a line, which is redshifted to near-IR 
wavelengths today (see, Santos et al. 2002). Thus, our mod- 
eling based on CLFs described above will not contain such a 
separate component especially if the Pop-Ill formation hap- 
pens at z > 8. The 3.5 /im IRB intensity associated with mod- 
els in the top and the bottom curve is <~ 2.5 and 0.3 nW m" 2 
sr" 1 , respectively. The angular power spectrum of clustering 
peaks at a multipole of ~ 10 3 , corresponding to the peak of 
the linear clustering when projected at these high redshifts. If 
the IRB excess is due to this component, the the clustering 
excess should peak at multipoles around £ 2 <~ 10 3 ; this is sim- 
ply a reflection of the projection associated with the peak of 
the linear mass power spectrum. As discussed in Cooray et al. 
(2004; Bock et al. 2006), in the H-band, the IRB intensity is 
roughly a factor of 10 higher given that the Pop-Ill spectrum 
is sharply rising towards lower wavelengths as the intensity 
there is dominated by the Lyman-a emission associated with 
recombinations during the reionization of the Universe. The 
two models shown in Figure 8 left panel bracket two extreme 
regions in terms of the reionization history, the density of first 
galaxies, their clustering strength or bias, and the Lyman-a 
photon production. It is clear that one should focus on degree 
to arcminute angular scales for anisotropy studies in the IRB 
instead of sub-arcminute angular scales considered by Kash- 
linsky et al. (2005) at multipoles above 10 4 . As is clear from 
Figure 8 left panel, such imaging need not be very deep since 
degree-scale shallow images that allow removal of sources 
down to a magnitude of 21 in the L-band can be easily used to 
probe the range suggested by the two extreme models. 

To briefly summarize our discussion so far, the excess clus- 
tering power spectrum measured in Kashlinsky et al. (2005) is 
consistent with faint galaxies with magnitudes between 22.5 
and 26 and predominantly at redshifts between 1 and 3 as ex- 
trapolated from measured LFs out to z ~ 1 . These faint galax- 
ies are unresolved in the deep Spitzer images, but not a signif- 
icant contributor to the total intensity of the IRB. If we were 
to explain the missing IRB with a new population of sources, 
such as Pop-Ill stars, then the clustering spectrum can be ex- 
plained, but there is an overall increase in fluctuations at mul- 
tipoles of <~ 10 3 which is not currently probed by unresolved 
clustering. 

6.3. Wide-Field Surveys for IRB Anisotropy Measurements 

While degree to arcminute angular scale fluctuation mea- 
surements allow first-stars at z > 8 to be more easily iden- 



tified, extending fluctuation measurements to low multipoles 
is challenging with Spitzer IRAC images alone due to both 
the small field-of-view and to problems such as uncertainties 
in relative calibration associated with stitching independent 
fields together to make a larger map. The Bootes field consid- 
ered here comes from a large number of independent images, 
but we are able to make clustering measurements of resolved 
sources since the catalogs are less affected by issues related to 
sensitivity variations from field to field. Nevertheless, we are 
also studying the possibility to use the 6 deg. 2 central region 
of the Bootes mosaic to at least place limits on the amplitude 
of unresolved source clustering at 10 arcminute angular scales 
to a degree angular scales. These results will be presented in 
an upcoming paper (Cooray et al. 2006). 

As mentioned in above, to separate various possibilities, the 
need for large-area clustering measurements is clear from the 
left panel of Figure 8 (see, also Figure 2 of Bock et al. 2006). 
While at multipoles of 10 4 to 10 5 the model clustering of a 
z > 8 diffuse component overlaps with faint galaxy clustering, 
at multipoles of 10 3 the clustering is well separated. Thus, 
at degree angular scales, instead of deep images one can use 
shallow images down to a brighter magnitude limit to search 
for an excess component than the magnitude limit one has to 
image when probing the same diffuse component with images 
that allow clustering studies only at arcminute angular scales 
and below. In consideration of the two issues discussed so far, 
large-angular scale excess and the large difference between 
measured IRB intensity and the predicted intensity based on 
counts at lower wavelengths than studied by Spitzer, we sug- 
gest observations that attempt to resolve the nature of the IRB 
excess should concentrate on shallow, wide-field imaging at 
wavelengths below 3.6 /im. 

The Cosmic Infrared Background Explorer (CIBER; Bock 
et al. 2006) is designed directly to address the "missing 
source" problem in IRB that is significant between 1 fim and 
2 fim through shallow, but wide-field images. CIBER uses 
wide-field imagers with a field-of-view of 4 deg. 2 in both I- 
and H-bands, supported by spectrometers to determine the 
contribution from zodiacal light and establish the spectrum 
of IRB intensity from 0.7 to 1.8 /im. With the wide-field 
coverage, CIBER can measure clustering of IRB light from 
multipoles below 10 2 to 10 4 covering the linear to non-linear 
regime from a single image. Furthermore, CIBER will image 
the Bootes field so that existing source counts from Spitzer 
and other ground-based imaging data can be used to remove 
point sources. The combination of resolved and unresolved 
source clustering by extending studies such as the one pre- 
sented here should allow us to pin down the sources respon- 
sible for the missing IRB and directly address the "missing 
source" problem of the IRB. 

In addition to CIBER, at wavelengths above 2.2 fim, the 
North Ecliptic Pole wide-field survey with the IR Camera 
on Akari (Matsuhara et al. 2006) will provide an additional 
dataset to study large-scale fluctuations in the near-IR back- 
ground. The shallow survey useful for clustering measure- 
ments will span over 6.2 deg. 2 and will allow 

clustering measurements from the linear to non-linear 
scales similar to the angular power spectrum established here 
for 3.6 /im resolved sources in the Bootes field. Finally, it 
could be that after the depletion of cryogen necessary for 
longer wavelength instruments, Spitzer will continue to op- 
erate only the first two channels of IRAC. Such a scenario 
could be exploited for a very wide-field survey which in return 
could be used to further improve the clustering measurements 
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FIG. 8. — IRB fluctuations from a z > 8 diffuse component that traces the linear density field (dot-dashed lines), galaxy clustering (solid lines), and shot-noise 
from galaxies (dashed lines that scale as r with increasing multipole). The measurements are from Kashlinsky et al. (2005) in the L-band with fainter symbols 
showing the background fluctuations after rebinning the original measurements. We consider three model descriptions for clustering using the CLF models (see 
text for details). The corresponding counts are shown and compared to completeness corrected measured counts in the right panel. In the right panel, in addition 
to total counts (solid lines), we also show counts of central (dashed lines) and satellite (dot-dashed lines) galaxies. These predictions for number counts are 
obtained by simply integrating the model LFs over the volume with luminosity converted to an observed flux as a function of redshift. The black lines show 
our default model and the shot-noise component can be roughly described as due to point sources with L band magnitudes between 22.5 and 26. This model 
also reproduces clustering suggesting that unresolved clustering spectrum measured by Kashlinsky et al. (2005) can be described by galaxies in the magnitude 
range between 22.5 and 26. If we populate IR galaxies down to a halo mass of 10 s M su n , we find counts shown by the blue line, with a slope consistent with the 
faint-end slope suggested in Savage & Oliver (2005). These counts lead to an order of magnitude higher shot-noise than measured with the clustering spectrum 
of IR anisotropics and a larger clustering amplitude than measured. Finally, with green lines, we only consider counts at z, > 5 for the same model, but in this 
case, we fail to explain the measurements given that these faint galaxies are not strongly clustered even at high redshifts. The two dot-dashed lines show a model 
expectation for clustering in the diffuse IRB from unresolved sources at z > 8 prior to complete reionization (Cooray et al. 2004). The 3.5 fim IRB intensity 
associated with models in the top and the bottom curve is ~ 2.5 and 0.3 nW irT 2 sr~' , respectively. The angular power spectrum of clustering peaks at a multipole 
of ~ 1 3 , corresponding to the peak of the lineal' clustering when projected at redshifts greater than 8. 



of both resolved and unresolved components. While the 
deep imaging data with Spitzer IRAC of the GOODS HDF-N 
and CDF-S fields have now allowed clustering measurements 
down to very faint magnitudes, these measurements are re- 
stricted to the non-linear regime. It will be very useful to 
measure clustering for sources at magnitude levels of 21 and 
fainter so that background anisotropy measurements at degree 
angular scales can be properly combined with resolved source 
clustering at these faint magnitude levels to obtain a complete 
picture of the IRB anisotropies. 

7. SUMMARY 

We have presented a measurement of angular power spec- 
trum of the clustering of near-IR sources at 3.6 /im in Spitzer 
imaging data of the GOODS HDF-N, the GOODS CDF-S, 
and the NDWFS Bootes field in several source magnitude 
bins. We also measured the angular power spectrum of re- 
solved sources in the Bootes field at K$ and /-bands using 
ground-based IR imaging data. In the three bands, /, Kg, 
and L, we have detected the clustering of galaxies on top of 
the shot-noise power spectrum at multipoles between £ ^ 10 2 
and 10 5 . The angular power spectra range from the large, 
linear scales to small, non-linear scales of galaxy clustering, 
and show a clear departure from a power-law clustering spec- 
trum for L-band galaxies when clustering in the Bootes field 
is combined with GOODS. We consider a halo model to de- 
scribe clustering measurements and establish the halo occu- 
pation number parameters of IR bright galaxies at redshifts 
around unity. 

The typical halo mass scale at which two or more IR galax- 
ies with L-band magnitude between 17 and 19 are found in 



the same halo is between 9 x 10 11 M Q and 7 x 10 12 M at the 
68% confidence level; this is consistent with the previous halo 
mass estimates for galaxies at z ~ 1 from clustering studies in 
surveys such as DEEP2. We have also discussed our results in 
the context of a recent measurement related to the unresolved 
IR background (IRB) anisotropies based on Spitzer imaging 
data by Kashlinsky et al. (2005). While the unresolved IRB 
fluctuations were measured at sub-arcminute angular scales, 
we have argued that the nature of suggested excess cluster- 
ing can be best studied with shallow, wide-field images that 
can make measurements of clustering from a few degree to ar- 
cminute angular scales with resolved L-band sources removed 
down to a magnitude level of about 21 and /-band sources re- 
moved down to a magnitude limit of 19. An attempt at making 
a wide-field image of the near-IR sky is planned with the Cos- 
mic Infrared Background Explorer (CIBER; Bock et al. 2006) 
at I- and H-bands, Akari at K-band and above, and potentially 
with Spitzer in the L-band. 
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